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Abstract: We investigate the role of asymmetries in the line spread function of the 2-degree field 
spectrograph and the variations in these asymmetries with the CCD, the plate, the time of observation 
and the fibre. A data-reduction pipeline is developed that takes these deformations into account 
for the calibration and cross-correlation of the spectra. We show that, using the emission lines of 
calibration lamp observations, we can ht the line spread function with the sum of two Gaussian functions 
representing the theoretical signal and a perturbation of the system. This model is then used to 
calibrate the spectra and generate templates by downgrading high resolution spectra. Thus, we can 
cross-correlate the observed spectra with templates degraded in the same way. 

Our reduction pipeline is tested on real observations and provides a significant improvement in the 
accuracy of the radial velocities obtained. In particular, the systematic errors that were as high as 
~ 20kms“^ when applying the AAO reduction package 2dfdr are now reduced to ~ 5kms“^. 

Even though the 2-degree Field spectrograph is to be decommissioned at the end of 2005, the analysis 
of archival data and previous studies could be improved by the reduction procedure we propose here. 
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1 Introduction 

The 2-degree held multi-object spectrograph at the 
Anglo-Australian Telescope can observe up to 400 ob¬ 
jects within a 2 degree held on the sky thro ugh 2 diher- 
ent sp ectrograph settings at the same time (iLewis et alJ 
l2002ll ^. Each spectrograph accepts light from 200 h- 
bres, which are positioned by a robot on a plate in the 
focal plane of the telescope, to produce spectra of low 
to medium resolution depending on the chosen grating. 
Two full sets of hbres on separate plates ensures that 
re-conhguring can be done in parallel with observing, 
hence minimizing time lost to the positioning of hbres 
during the night. 

With its high resolution gratings, 2dF should be 
an ideal instrument for studies of the kinematics of 
resolved stars, however, it has been shown that ra¬ 
dial velocities obtained using the 2dF spectrograph 
are plague d by systematic errors of the order of 10 — 
20 kms“^ Ic anno ^ I 2 OO 2 I : IStanford fc Cami^ 1200211 . 
While known, these systematic errors are not easily 
corrected for they are not constant over all the ob¬ 
servations. Observations performed on different CCDs 
or plates, during different nights or at different times 
during the same night can show variations in radial 
velocity that are worse than the 5 kms“^ internal pre¬ 
cision the 2dF should achieve. 

The most probable culprit for these variations is 
a change in the point spread function of the system 

^see also http://www.aao.gov.au/2df/ 
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across the CCDs due the optics of the cameras, tilts 
in the slit blocks or changes in the focus of the spec¬ 
trographs with temperature or orientation of the tele¬ 
scope, given that they are mou nted on the to p of the 
AAT close to the dome aperture llCa.nn onl2002fl . These 
deformations produce an asymmetric Line Spread Func¬ 
tion (LSF) on the extracted spectra which can be re¬ 
vealed by artihcially asymmetric emission/absorption 
lines. 

Our group recently undertook a radial velocity sur¬ 
vey of the region around the Canis Major (CMa) dwarf 
galaxy and other nearby overdensities in Red Giant 
Bran c h star s (RGBs) and Red Clump stars llMartin et alJ 
120041 . 12 OO. 5 II . Given its high number of targets on a 
wide held of view, the 2dF spectrograph is currently 
the best instrument for such a study. However, one of 
the key features of accretion streams is their very low 
dispersion in radial velocities. Thus, it is crucial that 
we achieve the lowest possible systematic errors on the 
velocities of our sample stars to be as close as possible 
to the internal dispersion of a population. Moreover, 
for the case of the Canis Major study, the velocity of 
the constituent star s is close to those of the thin disc 
llMartin et alJl^^4^ . In the regions where the CMa 
population is signihcantly contaminated by disc stars, 
a systematic error of 10 — 20 kms“^ could diffuse the 
signal of the extra-galactic population and could pre¬ 
vent detection. 

For these reasons, we constructed a reduction pipeline 
for 2dF observations that takes into account the asym- 
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metry of the LSF of the extracted spectra and cor¬ 
rects it. Since the shape of the LSF plays a role in 
the calibration of the observed spectra and their cross¬ 
correlation with templates to obtain radial velocities, 
our pipeline works as follows: 

1. determination of the deformations of the LSF 
across the CCD using the emission lines of the 
observed calibration lamps; 

2. use of the determined model of the LSF to cali¬ 
brate the arc spectra; 

3. use of the determined model of the LSF to gener¬ 
ate template spectra (from high resolution spec¬ 
tra) deformed in the same way as the observed 
spectra; 

4. cross-correlation of the observed spectra with 
the obtained template spectra. 

In section 2, we describe the different steps of the 
pipeline, while in section 3 the derived radial veloc¬ 
ities are compared to reference observations and we 
analyze the accuracy of the new velocities. Section 4 
summarizes the pipeline and the results. 

In the following, the 2dF observations are first cor¬ 
rected for the flat field and extracted using the 2dF 
Data Reduction pa ckage (2dfdr) provided by the AAO 
llTavlor et al.ll993l ^. Sky subtraction is also performed 
using the 2dfdr package. 

2 Observations and Reduc¬ 
tion 

The observations on which we apply our reduction 
pipeline in this section are part of our Canis Major 
radial velocity survey. We targeted the Red Giant 
Branch stars of a field centred on I — 240°, b = —8.8°. 

We employed two different spectrograph settings, with 
the 1200V grating on spectrograph 1 (covering 4600- 

O O 

5600 A at 1 A/pixel, these are the CCDl data) and 
with the 1200R grating on spectrograph 2 (covering 

8000-9000 A, also at 1 A/pixel, these are the CCD2 
dat a). 

IStanford fo GannonI ll2002li showed that observa¬ 
tions made with CGD2 are better than those obtained 
through CCDl. We will therefore concentrate of cor¬ 
recting the deformations of the LSF on CCDl. This 
is also justified by our analysis of CCD2 observations 
(see below, subsection 2.5). 

2.1 The LSF model 

For a given optical system, the output signal So(A) that 
is observed at wavelength A can be expressed as: 

s„(A) = (si*G)(A) (1) 

where Si is the input of the system and G is the line 
spread function of the system. When observed on a 

O 

CCD, each pixel X, of width S A receives part of this 
output signal, and has a value F{X) such that: 

^See also http://www.aao.gov.aU/2df/software.html#2dfdr 



Figure 1: 5105.541 A Cu emission line extracted 
from fibre 1 of a CCDl calibration lamp observa¬ 
tion. The emission line shows an asymmetry to¬ 
wards higher wavelengths that accounts for ^ 30 % 
of the total signal. 


j-X+O.SS 

F{X)= / So(A)dA. 

J x-o.ss 


( 2 ) 


When observing an Emission Line (EL) centred on 
fii, the input signal is: 


Si{X) = 6{X - fii) 
which produces the observed signal: 


(3) 


,(A) = 




G{r)S{X — fii — t) dr = G(A — /ii) . (4) 


This means that the signal, F{X), produced by this 
EL on pixel X of the 2dF with grating 1200V — for 

O 

which S = 1 A — is: 


^X-KO.5 


F{X) = 


G(A-Mi)dA. 


(5) 


If the system through which the observations are 
performed was perfect, the LSF of the system (G) 
would be a symmetric Gaussian function with a disper¬ 
sion that depends on the quality of the system. How¬ 
ever, this is not the case for the 2dF spectrograph. 
As can be seen on Figure 1 for an EL produced by 
the observation of a Copper-Argon calibration lamp 
on CCDl, the LSF is deformed and generates an EL 
with an asymmetric wing at higher wavelengths. This 
deformation is significant and can account for as much 
as 30% of the total signal of the EL. 

To model the deformed LSE, we represent it as 
the sum of two Gaussian functions. The first one, Gi, 
is the signal that would be observed through a perfect 
system and is only determined by the dispersion of the 
would be perfect system: cri. The second Gaussian 
function, G 2 , is considered as a perturbation of the 
system that produces the deformation. It is shifted 
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by Ap compared to Gi, has its own dispersion <72 and 
contains a fraction A of the total signal. Hence this 
model of the LSF can be expressed by: 


G(A) = (1-A)Gi(A)+AG2(A) 
1-A A 

O 2 \C71 J 




+ 


(71 




02 


(. 6 ) 


With this model and using equation |21 the value 
F{X) of pixel X of an emission line becomes: 


F{X) 


1-A 

pX+0.5 

\/^Ol j 

1 ^ 

X-0.5 

A 

pX-\-0.5 

02 

J X-0.5 


-1 / dA (7) 

_ 1 ^ - Am) ^ ^ 


dA. 


2.2 Fitting the model 

As expressed in equation Q the shape of an EL at a 
given position on the CCD is only determined by 5 free 
parameters: 

• pi: the theoretical centre of the emission line; 

• oi: the standard deviation of the theoretical, 
symmetric LSF; 

• A: the fraction of the total signal that is in the 
perturbation; 

• 02 '. the standard deviation of the perturbation; 

• Ap: the shift between the two Gaussian func¬ 
tions Gi and G 2 . Since the asymmetry of the 
emission lines is to higher wavelengths, Ap > 0. 


the centre of the region we will use to cross-correlate 

O O 

the spectra and the templates (the 4800 A to 5250 A 
region)^. We could determine the 5 free parameters 
of the model on the F function of each fibre, however, 
the low number of pixels over which F{X) ^ 0.0 (less 
than 20, see Figure 1) would produce high uncertain¬ 
ties on the parameters. Therefore, we choose to fit the 
same LSF model on up to 10 consecutive hbres (de¬ 
pending on the presence of dead fibres) at the same 
time. Indeed, the spectrograph is constructed with h- 
bres coming in blocks of 10 which means the hbres in 
each of these blocks should have analogous LSF® but 
changes may appear between each block. 

By extracting the EL of all the hbres and concate¬ 
nating them, we generate a new function F' that is 
modeled by F': 


« pX+O.S 

F'{X)^Y. G{X-pi)A\ ( 8 ) 

Jx-o.h 

In this new model, G is still the LSF model as ex¬ 
pressed in equation El Mi is the centre of the i*** EL 
and N is the number of hbres in the group (i.e. 10 
minus the number of dead hbres). This concatenated 
model F' now contains 4 + N free parameters: the N 
centres of the ELs {pi)i<i<N and the 4 parameters of 
the LSF {A, 01, 02 and Ap). This time however, there 
are 20 x A points to adjust these 4 -|- A free parame¬ 
ters which ensures lower uncertainties on the estimate 
of the parameters. 

Each group of hbres (1-10, 11-20, ...) is processed 
as follows: 

(i) for each hbre: 

• check if hbre is not dead; 


pi is the only free parameter that depends on the 
EL while the 4 others {A, 01 , 02 and Ap) dehne the 
shape of the LSE for the considered hbre, time and 
plate. By htting this model F of the EL on the EL 
observed on the CCD, F® (such as in Figure 1), we 
can deduce the parameters that dehne the LSF of the 
system for the hbre, wavelength, time and plate of the 
observation. 

To completely characterize the deformations of the 
LSF, it would be necessary to dehne these parameters 
in both the X (wavelength) and Y (hbre number) direc¬ 
tions of the CCD. However, since for our observations 
the region we use to cross-correlate the spectra and the 
templates is near the centre of the wavelength range of 
the CCD and since there are few ELs in this region, we 
only model the deformations of the LSE in the Y di¬ 
rection of the CCD. Yet, observations performed with 
different settings may require a complete two dimen¬ 
sion model of the LSF. 

Among the ~ 10 ELs that are observed with the 
1200V grating with a Copper-Argon calibration lamp, 
the best choice to ht the model F is the Cu line at 

O 

5105.541 A since it is the strongest line and it is at 

®In the following, the ~ functions (e.g. F) represent the 
observed data while the functions without the ~ represent 
the model to be fitted to the data (e.g. F). 


• extract the Cu line at 5105.541 A (around pixel 
500) to generate F; 

(ii) concatenate all the fibres of the group to gen¬ 
erate F'; 

(hi) fit model F' (4 + N parameters) to data F'. 

When extracting the EL, we take care to extract 
approximately the same region around the line for all 
the fibres, even if their centre varies with the posi¬ 
tion across the CCD. We detect the pixel Ac with the 
highest value of the EL and extract 6 pixels before Ac 
and 13 pixels after Ac to account for the asymmetry 
to higher wavelength in the shape of the EL. We then 
correct from a possible continuum by subtracting a lin¬ 
ear fit between the average value of the first 3 pixels of 
F and the average value of the last 3 pixels. Finally, 
to correct differences in the response of the fibres, we 
normalize each EL so that: 


^We choose this region for the cross-correlation with the 
templates because it is free of strong sky features. It is 
also located at the centre of the CCD along the wavelength 
axis, which should diminish potential asymmetries of the 
LSF generated along this axis. 

®This was tested by taking the fibres in smaller groups. 
The results are similar, but with higher uncertainties on 
the values of the parameters. 
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Pixel 

Figure 2: Concatenation of emission lines from a group of fibres (hollow circles, corresponding to the F' 
function). Since there are no dead fibres in this group, 10 emission lines were extracted and were used for 
the ht following the procedure described in the text. The best ht (the F' function) is shown as a thick 
line and faithfully represents the data. In particular, it reproduces the asymmetric, high wavelength wing 
of each emission line. 
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Figure 3: Evolution of the 4 LSF-related parame¬ 
ters A (bottom left panel), cti (bottom right), (T 2 
(top right) and Ap (top left) across CCDl. The 
three parameters representing the perturbation of 
the LSF {A, a 2 and Ap) all show substantial vari¬ 
ations across the CCD. The number of fibres, N, 
that were used for the fit is plotted as a thin dashed 
histogram in the bottom left panel. 


F{X)dX = 1.0 


(9) 


An example of the final data F' we obtain after the 
concatenation of the fibres of a group is shown as hol¬ 
low circles on Figure 2. 

To fit the model F' t o the data F', we use the 
mrqmin routine described in iPress et alJ lll992f> and de¬ 
termine the 4 + N parameters. Since determining all 
of the parameters at the same time proves difficult, we 
determine independently the N EL-related parameters 
{pi)i<i<N and the 4 LSF-related parameters A, ai, <72 
and Ap. The fitting procedure is as follows: 

(1) the N Pi values are determined for a loose grid 
of fixed values of A, ai, <72 and Ap. The {pi)i<i<N 
producing the lowest (in the sense of mrqmin) are 

used as input values; 

(ii) the {pi)i<i<N are fixed and the 4 LSF-related 
parameters are determined; 

(iii) these LSF-related parameters are fixed to their 
best values and used as input values to re-determine 
the {pi)i<i<N; 

(iv) these are once again fixed to re-determine the 
best values of A, ai, <72 and Ap. 

Steps (iii) and (iv) are used as a sanity check to 
verify the independence of the two sets of parame¬ 
ters. Since is reduced by less than 1% between 
step (ii) and step (iv), our hypothesis can be taken as 
valid and the 2 sets of parameters determined inde¬ 
pendently. Figure 2 shows the fit resulting from this 
procedure as a thick line. It is visibly good with a well- 
determined position for the ELs (the different pi) and 
it reproduces the asymmetric shape of the ELs with a 
large wing at higher pixel number. 


For a given calibration lamp observation, this pro¬ 
cedure is repeated for the 20 groups of up to 10 valid fi¬ 
bres. An example of the evolution of the 4 LSF-related 
parameters is shown on Figure 3. The first striking fea¬ 
ture is that the 4 parameters evolve across the CCD, 
which means that the asymmetry of the LSF is not 
constant for the different fibres of the same observation 
run. The most important parameter. A, which repre¬ 
sents the weight of the perturbation goes up to ~ 30% 
at the edges of the CCD but is as low as ~ 15% at it s 
centre. If this behaviour is expected llCa.nno 
the deformations at the edges of the CCD are so impor¬ 
tant they can easily explain strong systematic errors 
in the radial velocities if they are not accounted for. 
The parameters all have a roughly symmetric evolu¬ 
tion with values at the two edges of the CCD being in 
the same range and different from those in the centre. 

One of the assumptions in our LSF model is that 
it can be divided into a “theoretical signal” Gi (de¬ 
termined by (71) and a “perturbation” G 2 (determined 
by A, Ap and ( 72 ) that creates the asymmetry. This 
seems to be confirmed by the almost constant value of 
(71 that only evolves between 0.81 and 0.85 over the 
CCD, while at the same time, the “perturbation dis¬ 
persion”, (72, is in a range that is one order of magni¬ 
tude higher (between 1.3 and 1.7). Moreover, the value 
of (71 is similar to what is found on the symmetric ELs 
observed on CCD2 (see below, subsection 2.5). 

Though the parameters evolve in the expected way, 
a more detailed look at their behaviour shows many 
gaps and spikes where one would expect a smoother 
evolution over the CCD. This is partly due to changes 
in the number of fibres used to fit the model {N, dashed 
histogram in the bottom left panel of Figure 3): for in¬ 
stance, there are four dead fibres between fibre 81 and 
fibre 90 so, in the corresponding group, there remains 
only 6 ELs for the fit. For this group, the parame¬ 
ters Ap and (72 show a gap and have wider error bars. 
These gaps and spikes are also explained by the shal¬ 
low space which provides many slightly different 
LSF models with small differences in x^- 

Another interesting point to look at is the evolu¬ 
tion of the asymmetry on the LSF during a night of 
observation and the influence of the plate on which the 
observations are done. Figure 4 shows the evolution 
of the 4 parameters of the LSF for the 7 observation 
runs we made on 9th April 2004. For the same plate, 
the results are similar, with parameters that are really 
identical (e.g. (t\ and A) but the parameters of the 
“perturbation” tend to vary, even though they show a 
similar behaviour. Likewise, even if the deformation 
of the LSF evolves in a similar way on the two plates, 
there are some intrinsic differences between the two: 
for instance, the amplitude of the deformation tends 
to be higher at the center of the CCD on plate 1, the 
standard deviation of the “perturbation”, ( 72 , tends to 
be offset to lower values. 

Therefore, even if it is reassuring to observe a sim¬ 
ilar behaviour in the evolution of the deformation of 
the LSF during the night, it seems preferable to deter¬ 
mine a LSF model on the calibration lamp of each 
observation run. Using a global LSF model for all 
the observations done on a single plate during a single 
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Figure 4: Same as Figure 3 but with the evolution of the LSF parameters for all the observations that 
were made during a single night. The parameters for observations on plate 0 are shown on the left panels 
and those performed on plate 1 are shown on the right panels. Altough the parameters evolve in a similar 
way for all the fits, there are variations with plate and/or observation that prevent from using a generic 
set of parameters for the whole night. 


night nevertheless yields acceptable final results, even 
though the uncertainty of the derived radial velocities 
tend to be 20-30 percent higher than when using one 
LSF model for each calibration lamp. 

2.3 Calibration 

One part of the reduction where the asymmetry of the 
LSF has an important effect is the calibration of the 
spectra. Indeed, in a routine like identify in IRAF 
the centre, fio, of an emission line F defined over the 
range [Amin, Amax] is the wavelength at which the EL 
can be divided in two equally luminous parts. That is: 

_ p^max ^ 

/ F(A)dA= / F(A)dA. (10) 

Hence, the presence of the asymmetry in the EL tends 
to switch the “IRAF centre” fio to higher values than 
the “theoretical centre” fii which should be considered 
as a more valid position of the EL (see Figure 5). With 
the evolution of this offset directly related to the LSF 
parameter A, the difference /tq —changes with the fi¬ 
bre. And since we cross-correlate the observed spectra 
with templates that we generate (see next subsection) 
and do not have this calibration issue, the deduced 
radial velocity can have a systematic error of up to 
~ 10 kms"F 

To correct from this effect, we have to calibrate 
each arc spectrum using the LSF shape of the corre¬ 
sponding fibre, time and plate to determine the posi¬ 
tion of the “theoretical” EL (/ri) instead of the centre 
of the asymmetric, observed, EL (/ro). For a given fibre 
of a calibration lamp, we proceed as follows: 

(i) for each EL of the corresponding fibre on the 
calibration lamp, determination of its “theoretical cen¬ 
tre” fii by fitting the EL model of that fibre (equa- 
tion[71l with fii as the only free parameter (an example 


is given Figure 5); 

(ii) fit all the “theoretical centres” with a Legendre 
polynomial (~ 10 reliable ELs are used for this fit) to 
determine the wavelength of each pixel of the observed 
spectrum; 

(iii) the spectrum is re-binned so that each pixel is 

o o o 

1 A wide, starting at 4600 A and ending at 5600 A. 

Using the “theoretical centres” to calibrate the spec¬ 
trum ensures a calibration that is less altered by the 
shape of the 2dF LSF. Indeed, though we do not take 
into account possible evolutions of the LSF asymmetry 
along the wavelength direction of the CCD — we fit 
the same LSF model to all ELs — this technique gives 
a much better calibration than when using the sim¬ 
pler IRAF routine. To test that assumption one can 
measure the difference in radial velocity for the differ¬ 
ent fibres of a twilight observation. Since in such an 
observation, the spectrum on each fibre is that of the 
Sun, cross-correlating each calibrated spectrum with a 
reference spectrum (say, the one of fibre 100) directly 
reveals calibration-related systematic errors. This test 
was done for twilight observations of 9th April 2004 
and the difference in radial velocity between the differ¬ 
ent fibres is shown on Figure 6 for a calibration using 
the IRAF identify routine (left) and the procedure 
described above (right). 

It is clear that directly using the IRAF routine 
is not satisfactory with a velocity difference between 
fibre 100 and the other fibres that is not constant and 
which follows an arc with increasing discrepancies at 
the edge of the CCD. On the other hand, this effect 
is corrected when using our calibration procedure and 
all the fibres produce a similar radial velocity, with a 
small dispersion of 2.8kms“^. 
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Pixe 

Figure 5: Difference in the central position of an 
EL for two calibration methods, po is defined as 
in equation nni and corresponds to the value given 
by the IRAF routine identify while pi is the cen¬ 
tre of the “theoretical signal”, Gi, of the LSF, as 
defined in equation El The observed values, F, of 
the EL are shown as hollow circles, while the LSF 
model, F, determined following the procedure of 
subsection 2.2 is represented as a thick line. Its 
two components (the “theoretical signal” Gi and 
the “perturbation” G 2 ) are plotted as thin dashed 
lines. In this example, there is a shift of 0.28 pixels 
between the centres of the LSF, po and /ri, deter¬ 
mined by the two different methods. 


2.4 Generating templates 

Measuring the radial velocity of a star requires one to 
cross-correlate the stellar spectrum with that of a tem¬ 
plate star of known radial velocity. Generally, a tem¬ 
plate spectrum is obtained by observing, along with 
the target stars, a few well-known stars. However, due 
to time constraints, one is forced to observe a small 
number of templates during a specific night, with a 
specific fibre of the 2dF, whereas as we saw earlier, 
the asymmetry of the LSF can evolve with time and 
certainly changes with the fibre (see Figure 4). Con¬ 
sequently, the cross-correlation would be produced by 
two spectra deformed in different ways, which can pro¬ 
duce systematic errors on the measured shift between 
the two spectra of order of 0.2 pixels (corresponding to 
as much as 15kms“^ on CCDl). Hence, it is necessary 
to take the LSF shape of the templates into account 
when doing the cross-correlation. 

Since it is cumbersome to deconvolve the spectrum 
of the templates and then reconvolve it with the LSF 
of the fibre one is interested in, we choose to gener¬ 
ate artificial templates by convolving High Resolution 
(HR) spectra by the deformed LSF. This way, the tar¬ 
get and template spectra are similarly deformed and 
their cross-correlation no longer has systematic effects 
due to differences in the deformation of the LSF of the 
two spectra. For each group of 10 fibres, we proceed 
as follows: 

(i) convolution of the HR template with the corre- 


o IRAF routine Section 2.3 calibration 



-10 -5 0 5 10 -10 -5 0 5 10 
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Figure 6: Difference in radial velocity between the 
different fibres and fibre 100 for a twilight obser¬ 
vation. On the left panel, the spectra were cal¬ 
ibrated using the IRAF routine identify while 
on the right panel, the calibration was performed 
using the shape of the LSF, as explained in subsec¬ 
tion 2.3. The IRAF calibration generates an arc 
shaped systematic error that depends on the fibre 
number while our procedure does not introduce 
such errors. 


spending LSF (as defined in equationEland determined 
as explained in subsection 2.2); 

(ii) downgrading the HR, convolved template to 

O 

a low resolution, 2dF-like template, (1 A/pix for the 
1200V grating); 

(iii) cross-correlation of the observed spectra with 
the artificial template using the f xcor routine in IRAF. 

We choose the HR spectra from the UV ES Paranal 
Observatory Project^ iBaenulo et al.lf2003li which has 
the advantage of providing the spectrum of stars of 
numerous spectral types (from type O stars to type M 
stars) and hence provides the opportunity to choose 
a template that resembles the observed stars. More¬ 
over, the UVES observations made through the 580L 

O 

grating (4760-5770 A) correspond to the range we are 
interested in (4800-5250 A). 

2.5 CCD2 

As is visible on Figure 7, the calibration lamps ob¬ 
served on CCD2 show emission lines that do not have 
asymmetries like CCDl. In fact, the emission lines 
are well fitted by a single Gaussian model (the thick 
line of Figure 7) with uncertainties on the position and 
dispersion of the Gaussian lower than 0.01 pixel. It is 
therefore not necessary to apply the previous reduc¬ 
tion techniques to obtain precise radial velocities from 
observations on CCD2. However, the symmetry of the 
emission lines does not mean that there are no issues 
with this CCD. Indeed, analyzing the same EL for all 

the fibres (the 8408.2096 A CuAr line roughly located 

^ http://www.sc.eso.org/santiago/uvespop/index.html 
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Figure 7: 8408.2096 A CuAr emission line ex¬ 
tracted from fibre 1 of a CCD2 calibration lamp 
observation (hollow circles). The emission line is 
symmetric and can be fitted by a single Gaus¬ 
sian function of internal dispersion cr = 0.74 (thick 
line). 


Figure 8: Evolution of the internal dispersion, cr, 
of the LSF across CCD2 when fitted by a single 
Gaussian function as in Figure 7. Although the 
LSF remains symmetric across CCD2 (contrary to 
GCDl), its width varies substantially. 


at the centre of the CCD) shows there are substantial 
changes in the internal dispersion, cr, of the Gaussian 
fit. Once again, a has an arc-like evolution across the 
CCD (see Figure 8). Moreover, the evolution differs 
with the plate on which the observations were per¬ 
formed. Even though an increase in the dispersion is 
expected at the edges of the CCD, the magnitude of 
the variation (from a = 0.6 at the centre up to cr = 0.75 
at the edge of plate 0 observations or cr = 0.85 at the 
edge of plate 1 observations) is troublesome. In addi¬ 
tion, there is at least one group of fibres (the 90-100 
group on plate 0) that has a higher internal dispersion 
than its neighbouring groups. 

Altough the symmetry of the ELs means the 2dfdr 
reduction package should suffice to reduce the data 
precisely, the variation of the internal dispersion of the 
LSF on CCD2 is to be kept in mind if some systematic 
effects remain at the end of the reduction. 

3 Results 

To analyze the efficiency of the procedure we present 
here, we observed one of the fields of our Canis Major 
survey twice. The hrst observation run was performed 
on 9th April 2004, using the spectrograph 1 settings 
(CCDl and 1200V grating; this is our April dataset) 
and the second run was performed on 3rd December 
2004, this time using spectrograph 2 settings (CCD2 
and 1200R grating; this is our December dataset). As 
we showed, the radial velocities derived from the De¬ 
cember observations on CCD2 should not have impor¬ 
tant systematic effects. Moreover, the use of the 1200R 
grating produces low uncertainties on the velocities 
(~ 2kms“^). In the following, we will use the De¬ 
cember radial velocities as a reference and analyze the 
influence of our pipeline on the April radial velocities. 

When using the 2dfdr package to reduce the April 


observations (left panels of Figure 9), the radial ve¬ 
locities are broadly ~ 10 — 20 kms“^ lower than the 
reference velocities measured in the December run, but 
some stars have offsets as low as —30 kms“^ or as high 
as 20 kms“^. At the same time, the uncertainties on 
these velocities are mainly in the range 6 — 8 kms“^ 
with, here as well, groups of outliers with uncertainties 
increasing to ~ 15 kms“^ on different portions of the 
CCD (for low hbres or hbres between 130 and 150). 

In contrast, the reduction of the April observations 
on CCDl using the pipeline described in this paper 
(right panels of Figure 9) yields lower offsets compared 
to the reference velocities. The mean offset for all the 
CCD has changed from —9± 10kms“^ to 0±8kms“^ 
and there are only a few outliers (~ 2). Similarly, the 
uncertainties on each radial velocity are slightly lower 
and are more clustered, with no more groups of hbres 
producing outliers with high uncertainties. 

Though the use of the reduction pipeline described 
in this paper produces much more reliable radial veloc¬ 
ities, low systematic offsets still remain. On average, 
this effect is limited (~ 5 kms“^ and certainly less 
than 10 kms“^) and could either be due to the reduc¬ 
tion pipeline that may not completely correct the LSF 
asymmetry effects (but it has been shown that there 
are no more calibration-related issues) or it could in 
fact be due to systematic errors in the reference ve¬ 
locities. We have shown that the emission lines ob¬ 
served through CCD2 do not present any asymmetry 
but what remains unclear is the role that the varia¬ 
tions of the width of the LSF of CCD2 (see Figure 8) 
has on the radial velocities. 


4 Summary 

We have constructed an improved reduction pipeline 
for 2dF observations. The systematic offsets that ap¬ 
peared between radial velocities measured from CCDl 
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Figure 9: Difference in radial velocities between the April and December observations when reducing 
the data using the 2dfdr package (left panels) and the reduction pipeline presented in this paper (right 
panels). For each set of panels, the offset in radial velocity, Avr = fr.Apr ~ i^r.Dec, has been plotted (left) 
along with the uncertainties on these measurements (right). Using the reduction techniques developed in 
this paper produces lower systematic errors with Avr 0 kms“^ and a lower internal dispersion of the 
offsets. 


data and those measured from CCD2 data can now be 
corrected. It has been shown that they are due to an 
asymmetric line spread function, which for the data we 
analysed, particularly affects CCDl. This effect can be 
countered by: 

1. modeling the asymmetric LSF by the sum of two 
Gaussian functions — the first one represents 
the LSF of the system if it were perfect and the 
second one is a perturbation of this system; 

2. fitting the LSF model to emission lines of the 
calibration lamps for each group of 10 fibres; 

3. using this LSF model to calibrate the correspond¬ 
ing spectra; 

4. using this LSF model to downgrade high reso¬ 
lution template spectra so that the observations 
can be cross-correlated with templates similarly 
deformed. 

The difference in radial velocity in a set of Red Gi¬ 
ant Branch stars that has been observed twice (once 
on each CGD) is highly reduced when using this pro¬ 
cedure (~ ±5 kms“^). Moreover the radial veloc¬ 
ity measurements show a more consistent behaviour 
throughout the CGD than when using only the AAO 
reduction pipeline. 

Even though the 2dF spectrograph is to be decom¬ 
missioned at the end of this year, archival data from 
the past decade could be improved by the reduction 
procedure we propose here. In particular, studies of 
the kinematics of Galactic populations could benefit 
highly from a reduction of the internal dispersion of the 
radial velocity datasets. The pipeline has only been 
shown to work for a particular dataset, yet we suspect 
that observations obtained through other gratings or 
with other central wevelengths and that present the 
same kind of systematic effects can be corrected in a 
similar way. 
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